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ABSTRACT 



Heat transfer performance of horizontal, integral-fin tubes made of copper, 
aluminum, copper-nickel, and stainless steel was evaluated using a boiler and steam 
condenser assembly. Testing was done at vacuum and atmospheric pressure conditions. 
The tubes tested had an inner diameter of 12.7mm, a root diameter of 13.88mm, and fin 
heights ranging from 0.5mm to 1.5mm, in 0.25mm increments. The outside heat transfer 
coefficient was determined first by finding the overall heat transfer coefficient, Uq, then by 
using the Modified Wilson Plot Technique. 

The results indicated that the performance of a finned tube is very dependent on fin 
height and tube material. Moreover, the results were compared with the predictive models 
of Beatty and Katz, Rose, Adamek and Webb, and Honda et al., with a modified version 
of the Rose model demonstrating the best predictive capabilities. 
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NOMENCLATURE 



effective surface area as defined by eqn. (5) , 
surface area of fin flank as defined by eqn. (6) , 
surface area of fin tip as defined by eqn. (7) , itt^ 
inside surface area of test tube, in^ 
outside surface area of smooth tube m^ 
outside area of test tube for one pitch length, m^ 
unfinned surface area as defined in eqn. (8) , m^ 
constant used by Rose [Ref. 4], equal to 2.96 

constant used by Rose [Ref. 4], equal to 0.143 

constant used by Rose [Ref. 4], equal to 0.143 

constant used by Rose [Ref. 4], equal to 0.143 

assumed leading coefficient for hj as in eqn. (25) 
specific heat at constant pressure, J/(kg K) 
equivalent diameter as defined in eqn. (3) , m 
inside diameter of test tube, m 

outside diameter of test tube, or smooth tube, m 
root diameter of finned tube, m 

fraction of unflooded fin flank surface area that is 
covered with condensate 

fraction of unflooded interfin surface area that is 
covered with condensate 
gravitational constant, 9.81 m/s^ 
specific enthalpy of vaporization, J/kg 
inside heat transfer coefficient, W/ (m^ K) 
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outside heat transfer coefficient, W/ (m^ K) 

thermal conductivity, W/(m K) 

thermal conductivity of coolant, W/ (m K) 

thermal conductivity of condensate film, W/(m K) 

constant as defined in eqn. (28) 

constant as defined in eqn. (29) 

length of test tube, m 

fin flank length as defined in eqn. (4) , m 

log mean temperature difference, K 

mass flow rate of coolant, kg/s 

number of fins per unit length of tube, m*’ 

Prandtl number 

fin flank heat flux as defined in eqn. (10) , W/m^ 

interfin heat flux as defined in eqn. (11), W/m^ 

fin tip heat flux as defined in eqn. (12) , W/m^ 

heat transfer rate as defined in eqn. (19) , W 

Reynolds number 

interfin spacing, m 

fin thickness, m 

coolant inlet temperature, K 

coolant outlet temperature, K 

film temperature, K, or constant as in eqn. (16) 

steam temperature, K, or constant as in eqn. (17) 

steam saturation temperature, K 

constant as defined in eqn. (15) 

tube outside wall temperature (at fin base) , K 
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overall heat transfer coefficient, W/ (m^ K) 



U 



o 



GREEK SYMBOLS 



a assumed leading coefficient to find 

AT temperature difference across the condensate film, K 

fin efficiency 

6 constant as defined in eqn. (27) 

enhancement ratio for a given temperature difference 
as defined in eqn. (14) 
dynamic viscosity, kg/(m s) 

condensate film dynamic viscosity, kg/ (m s) 
p density, kg/m^ 

condensate film density, kg/m^ 
p^g fluid/vapor density difference, kg/m^ 

py vapor density, kg/m^ 

4) condensate flooding angle as defined in eqn. (13) 

o condensate surface tension, N/m 

5(4)) constant as used in eqn. (11) 

Q Petukhov-Popov function as defined in eqn. (26) 
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I. 



INTRODUCTION 



A. BACKGROUND 

Today, all over the world, steam plants are being used to 
provide power and electricity on land, and to propel ships and 
submarines at sea. Because of this extensive use of steam 
plants in general, and condensers in particular, it becomes 
apparent that any enhancement in the performance of a 
condenser could be of enormous benefit. For example, 
electricity could be generated cheaper, fuel consumption could 
be reduced, or ship speeds could be increased for a given 
power plant. 

One method of increasing condenser, and hence steam plant 
performance, is to use "enhanced" condenser tubes. These tubes 
offer an increase in performance by enhancing the heat 
transfer on either the inside or outside of the tubes. 
Therefore, using these tubes would allow for smaller, more 
efficient future condensers. Moreover, higher efficiency could 
be achieved for existing power plants by retubing with 
enhanced tubes. 

One type of enhanced tube is the integral-fin tube. An 
integral-fin tube is a tube with circumferential fins on its 
outside, manufactured by machining the material between the 
fins away. As the fin material always was part of the original 
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tube stock, there is no contact resistance between the fin and 
the tube wall, (ie, The fin is an integral part of the tube.) 

There are two main reasons why integral-fin tubes are 
enhanced over smooth tubes. One reason is because of the added 
surface area presented by the fins for heat transfer. The 
other reason is the interaction between the surface tension of 
the condensate and the fins themselves. 

Increasing the surface area of a tpbe, one might surmise, 
would be very important in enhancing the heat transfer 
performance of a tube. After all, the more surface area there 
is, the more area there is for heat transfer. However, one 
would also surmise that there must be a limit to heat transfer 
enhancement. Particularly with lower conductivity materials, 
it is intuitively obvious that there is a fin height beyond 
which no further practical heat transfer increase will occur. 
This limit in heat transfer rate results from the competitive 
effects of increased condensing surface, and decreased heat 
conduction (fin efficiency) through the fin as fin height 
increases. The effect of fin efficiency during single phase 
heat transfer is well known in setting a proper integral fin 
height. 

The interaction between the fins and the condensate 
during condensation is a complex one, with two competing 
effects arising from surface tension. One effect is to thin 
the condensate film on the upper part of the tube. This is 
called the unflooded region. On the lower part of the tube. 
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the presence of the fins causes condensate to be retained in 
the space between the fins. This is called the flooded region. 
These regions are shown in Figure 1. 

The unflooded region demonstrates enhanced heat transfer. 
This is because the condensate film on the tube wall and fin 
flanks is kept very thin by the action of surface tension and 
gravity. As the condensate has a much lower thermal 
conductivity than the typical metal tube, its thinning 
increases the amount of heat transfer. 

Again, because of the low conductivity of the condensate, 
the heat transfer is drastically reduced in the flooded 
portion of the tube. When compared to the unflooded portion, 
the amount of heat transfer provided by the flooded portion is 
very small. 

Unfortunately, by increasing the fin height, the flooded 
portion of the tube is increased as well, again because of the 
effects of surface tension. This tends to reduce the amount 
of heat transfer, and competes directly with the enhancing 
factor of increased tube surface area mentioned earlier. 

Much work has already been done with integral-fin tubes at 
the Naval Postgraduate School (NPS) and elsewhere. However, 
the vast majority of work has been done with copper tubes 
because of its high thermal conductivity and ease of 
fabrication. Because of strength and/or corrosion concerns. 
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Figure 1 Schemacic ot Condensate Retention Angle on 

Finned Tubes and Condensate Wedge (illustrated 
by the gray sections) 



most condensers use tubes made of copper-nickel, bronze, 
stainless steel, or titanium, all of which have much lower 
thermal conductivities than copper. 

B. PREDICTIVE MODELS 

It is obvious that enhanced tubes are advantageous. 
However, being able to predict their performance would be even 
more advantageous. After all, how does one design a condenser 
when the performance of the tubes isn't well known? For that 
matter, how does one tell if performance of enhanced tubes is 
worth the added cost of manufacturing them? 

Nusselt [Ref. 1], in 1916, was the first to successfully 
predict the performance of smooth tubes. Since then, Beatty 
and Katz [Ref. 2], Adamek and Webb [Ref. 3], Rose [Ref. 4], 
and Honda et al. [Ref. 5] have all attempted to predict, with 
varying degrees of success, the performance of integral-fin 
tubes . 

There is very little experimental validation of the 
previously mentioned integral-fin models and virtually all the 
data are with copper tubes (though Jaber and Webb [Ref. 6], 
have done some very recent work with other materials) . 
Therefore, the previously mentioned models remain essentially 
unproven with regard to tubes that would be used in actual 
condensers . 
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C. NAVAL POSTGRADUATE SCHOOL CONDENSATION RESEARCH 

This thesis is part of an ongoing research program to 
study enhanced condensation. Much work has been done over the 
years with integral-fin tubes of various dimensions, though 
most has been done only with copper tubes. Mitrou [Ref. 7], 
and most recently Cobb [Ref. 8], looked at tubes of different 
materials but with only limited variations of fin height. 



D. OBJECTIVES 

The main objectives of this thesis are as follows: 

1. Obtain repeatable data for integral-fin tubes made of 
different materials, to study the effects of thermal 
conductivity on tube performance. 

2. Compare data for tubes of the same material but different 
fin heights, to demonstrate the effect of fin height on tube 
performance. 

3. Compare the experimental results with available 
predictive models, to validate the models. 
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II. A REVIEW OF RELEVANT PREDICTIVE MODELS 



A. NUSSELT MODEL 

As mentioned previously, Nusselt [Ref. 1] was the first to 
formulate an equation for the average heat transfer 
coefficient for a smooth horizontal tube during film 
condensation : 



h^=0 .728 



f(pf-p y) 



( 1 ) 



In order to develop his equation, Nusselt assumed that 
the tube operates in a quiescent vapor, that is a vapor with 
zero velocity. While his model remains generally valid, in 
reality any vapor in a condenser will have some velocity. 
Assuming downward flow, the vapor velocity would tend to thin 
the condensate film and enhance the heat transfer above what 
the Nusselt model predicts. 

B. BEATTY AND KATZ MODEL 

In 1948, Beatty and Katz [Ref. 2] formulated an equation 
for the average heat transfer coefficient for integral-fin 
tubes. They took into account the thermal conductivity of the 
wall material in order to accurately model the effect of the 
fins. However, to simplify the problem they neglected the 
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effects of condensate surface tension. For rectangular shaped 
fins, their equation takes the form: 



where 



ho=0.689 






n/4 



( 2 ) 



n /4 



D 



=1 .3ti 



e< 7 j 



^fs 









KfD. 



X/4 



(3) 



and 



L = 7l ■ 



4£), 






(4) 

(5) 






rifV. (dI-dI) 



Af^=n^%D^t 



( 6 ) 

(7) 



A^=nfUD^s 



( 8 ) 



As Beatty and Katz ignored surface tension, one would 
expect their model to perform better for low surface tension 
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fluids, such as refrigerants, than it would for water. Also, 
the model would predict the performance better under, high 
pressures and hence, high saturation temperature conditions 
where surface tension would be lower. 

C. ROSE MODEL 

Rose [Ref. 4] in 1993, developed a simple but complete 
model for determining the outside heat transfer coefficient 
for integral-fin tubes. Unlike Beatty and Katz [Ref. 2], he 
took into account the effects of surface tension, gravity 
induced drainage from the tube, and condensate flooding. He 
did, however, choose to ignore the effects of fin efficiency 
as he primarily dwelled on copper tubes which have a very high 
fin efficiency. Rose's equation for the outside heat transfer 
coefficient for an integral-fin tube is; 



^o= 



TiD^tg^ 



71 









ATA 



( 9 ) 



tot.p 



where q^, q^, and q^ are the heat fluxes from the fin flanks, 
inter fin space, and fin tips: 



(Jr 





0.943‘Pf^S^ a] 




K 



1/4 



( 10 ) 
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{5{<t>) a] 




°r • s’ J ] 



1/4 



( 11 ) 



and: 



Qt= 





■0.724pf^g o] 




■Do 



1/4 



( 12 ) 



and the condensate flooding angle <|) is 

<j)=cos‘^ 



4o 






-1 



(13) 



The quantities and represent the fraction of the 
unflooded portion of the interfin space and the fin flanks 
that are flooded with condensate. 

Moreover, Rose defines the enhancement ratio as the 
ratio of the predicted outside heat transfer coefficient for 
a finned tube to that predicted by Nusselt at the same film 
temperature difference. This ratio is given as: 









Dj.{s+t) 






dI-dI) 



2Dj.(s+C) 



r „ SI (s+c) ® 



(14) 



It 



where : 









OD^ 



11/4 



0.728"Pf^grt\ 



( 15 ) 
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T - [ Q-943 ]^ Dr 
^ I 0.728] 0 .728V 



1/4 



( 16 ) 



and : 



0.728^ ^ 0 .728*p 



1/4 



T = 

S 



(17) 



Note that these equations contain four unknown 
coefficients, Bl, Bs, Bf, and Bt. Rose curve fitted these 
equations to existing experimental data for copper tubes at 
atmospheric pressure (only) and determined that Bl should be 
2.96, while Bf, Bs, and Bt, were all equal to 0.143. 

Cobb [Ref. 8], in 1993 modified the Rose model to include 
the effects of fin efficiency. The modified Rose model 
therefore takes the form of: 










coc,p 
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D. ADAMEK AND WEBB MODEL 



Adamek and Webb [Ref. 3] use a far different approach to 
determine the outside heat transfer coefficient. Like Rose 
[Ref. 4], gravity drainage, surface tension and the flooding 
angle are all taken into account. However, that is where the 
similarity ends. 

Adamek and Webb chose not to ignore the effects of fin 
efficiency. Furthermore they decided to look at a length of 
tube which stretches from the midpoint at the tip of a fin to 
the midpoint of its adjacent interfin space (see Figure 2) . 
The surface between those two points is then broken up 
into eight discrete segments, namely, ba, aO, 01, 12, 23, 34, 
45, and 56. For each of these segments, a local condensation 
rate for the condensate surface is calculated. These 
condensation rates are then summed for both the flooded and 
unflooded portions of the tube. In addition, condensate film 
thicknesses are determined for each of the eight segments. The 
outside heat transfer coefficient is then a function of the 
condensation rates, film thickness, fin efficiency, 
temperature difference, and enthalpy. A major disadvantage of 
this model is its complexity compared to the models of Rose 
[Ref. 4] or Beatty and Katz [Ref. 2], and a numerical solution 
is required to solve the problem. 
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Figure 2 Half Fin/Interfin Space as 
Analyzed by Adamek and Webb 



13 



E. HONDA MODEL 



The Honda et al. model [Ref. 5], like that of Adamek and 
Webb is quite complex, but is the most comprehensive model 
available. Like Adamek and Webb [Ref. 3], the condensate film 
thickness is calculated, and gravity and surface tension 
effects are considered. For Honda's model, three cases are 
considered based on fin spacing and condensation rate (see 
Figure 3 from Ref. 5) . Different sub-models are used for each 
case. These cases are a function of fin spacing and 
condensation rate and are used because it is expected that the 
depth of the condensate film in the inter-fin space would have 
a significant impact on the amount of heat transferred. 

Honda et al. [Ref. 5], however, take into consideration 
the properties of the test tube coolant, the inside heat 
transfer coefficient, and the tube wall conductivity in 
analyzing the heat transfer from the vapor to the coolant, and 
then determine the temperature field in the tube and fins. 
Therefore, their predicted outside heat transfer coefficient 
is a function of coolant properties, inside heat transfer 
coefficient, tube wall conductivity, fin efficiency, film 
thickness, and surface tension and temperature difference. 
This comprehensive analysis, however, requires a numerical 
solution. 
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ricjure 3 Three Sub-Cases of the Honda et al. Model 



Ill . EXPERIMENTAL APPARATUS 



A. SYSTEM AND SYSTEM INSTRUMENTATION OVERVIEW 

The system apparatus and instrumentation are identical to 
that as described by Cobb [Ref. 8]. A major computer upgrade 
is in progress, but has not yet been installed. 

B. TUBES TESTED 

As mentioned in the introduction, little experimental work 
has been done with tubes made of materials other than copper. 
For this work, tubes made of copper, aluminum, 90/10 copper- 
nickel, and 316 stainless steel were used in order to 
determine the relationship between tube heat transfer 
performance and tube thermal conductivity. The thermal 
conductivities for the tubes used were curve-fitted by Cobb 
[Ref. 8] for the temperature range of this work, from data 
taken from [Ref. 9]. Table I lists the thermal conductivities. 
TABLE I. THERMAL CONDUCTIVITIES OF TUBE MATERIALS 



1 MATERIAL 


THERMAL CONDUCTIVITY 
(w/ (m K) ) 


COPPER 


390.8 


ALUMINUM 


231.8 


COPPER-NICKEL 


55.3 
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THERMAL CONDUCTIVITY 
(w/(m K)) 


STAINLESS STEEL 


14.3 



All tubes tested contained a heatex insert. The heatex 
insert is an insert of wire loops and is used to promote 
repeatable, consistent, turbulent flow on the inside of the 
tubes to enhance the inside heat transfer coefficient and 
lower the inside thermal resistance. The tubes tested, and 
their dimensions are listed in Table II. 



TABLE II. SPECIFICATIONS FOR TUBES TESTED 



TUBE 

MATERIAL 


ROOT 

DIA. 

(MM) 


FIN 

HEIGHT 

(MM) 


OUTER 

DIA. 

(MM) 


FIN 

THICKNESS 

(MM) 


FIN 

SPACING 

(MM) 














COPPER 


13.88 


1.50 


15.88 


1.00 


1.50 


COPPER 


13.88 


1.25 


16.38 


1.00 


1.50 


COPPER 


13.88 


l.SO 


15.88 


1.00 


1.50 


COPPER 


13.88 


0.75 


15.38 


1.00 


1.50 


COPPER 


13.88 


0.50 


14.88 


1.00 


1.50 


COPPER 


13.88 


SMOOTH 


13.88 
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TUBE 

MATERIAL 


ROOT 

DIA. 

(MM) 


FIN 

HEIGHT 

(MM) 


OUTER 

DIA. 

(MM) 


FIN 

THICKNESS 

(MM) 


FIN 

SPACING 

(MM) 


ALUMINUM 


13.88 


1.50 


16.88 


1.00 


1.50 


ALUMINUM 


13.88 


1.25 


16.38 


1.00 


1.50 


ALUMINUM 


13.88 


1.00 


15.88 


1.00 


1.50 


ALUMINUM 


13.88 


0.75 


15.38 


1.00 


1.50 


ALUMINUM 


13.88 


0.50 


14.88 


1.00 


1.50 


ALUMINUM 


13.88 


SMOOTH 


13.88 






COPPER- 


13.88 


1.50 


16.88 


1.00 


1.50 


NICKEL 


\ 










COPPER- 


13.88 


1.00 


15.88 


1.00 


1.50 


NICKEL 












COPPER- 


13.88 


0.75 


15.38 


1.00 


1.50 


NICKEL 












COPPER- 


13.88 


0.50 


14.38 


1.00 


1.50 


NICKEL 












STAINLESS 


13.88 


1.50 


16.88 


1.00 


1.50 


STEEL 












STAINLESS 


13.88 


1.25 


16.38 


1.00 


1.50 


STEEL 
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TUBE 

MATERIAL 


ROOT 

DIA. 

(MM) 


FIN 

HEIGHT 

(MM) 


OUTER 

DIA. 

(MM) 


FIN 

THICKNESS 

(MM) 


FIN 

SPACING I 

(MM) 


STAINLESS 

STEEL 


13.88 


1.00 


15.38 


1.00 


1.50 


STAINLESS 

STEEL 


13.88 


0.75 


14.88 


1.00 


1.50 


STAINLESS 

STEEL 


13.88 


0.50 


14.38 


1.00 


1.50 



19 



IV. EXPERIMENTAL PROCEDURES AND DATA ANALYSIS 



A. SYSTEM OPERATION AND TUBE PREPARATION 

System (see Figure 4) operation was identical to that 
given by Cobb [Ref. 8]. For both atmospheric and vacuum runs, 
non-condensable gasses were removed by use of a vacuum pump. 
Simultaneously, the boiler heaters were turned on, and flow 
was initiated in the test tube. Once steady conditions were 
reached for the vacuum (saturation temperature of 48.7 degrees 
C) or atmospheric (saturation temperature of 100.0 degrees C) 
runs, cooling water flow was adjusted to 80% in the test 
tube. 

At this point data collection commenced. The data 
collection procedure was repeated and the temperatures checked 
for consistency before saving them. If the data were 
sufficiently consistent, (+/- 1%) the flow through the test 
tube was repeated with the flow meter reduced to 70%. This 
process continued down to 20% flow in the test tube and was 
then repeated from 20% back up to 80%. 

Tube preparation was also identical to that given by Cobb 
[Ref. 8] with the following exception: 

• For aluminum tubes only, the treatment was stopped once a 
continuous oxide layer has been formed on the surface of 
the tube, but before dimensional changes had occurred 
because of excessive corrosion due to the high reactivity 
of aluminum. 
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Figurs 4 Schematic of the Single Tube Tesu Apparatus 
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B. COMPUTER CODES 



Three different computer codes were used for analysis in 
this work. The first of these codes was used to take the raw 
data and do initial processing, while the second and third 
were codified versions of the previously mentioned predictive 
models . 

1. DRPALL 

"DRPALL" is the name of the data acquisition and 
initial processing program. It is an HPBASIC program and 
remains unchanged from that described by Cobb [Ref. 8]. 

When used, the DRPALL program asks the user for information 
regarding test tube material type and configuration. Once the 
operator is ready to commence data taking, DRPALL either 
measures directly via an HP 3497 Data Acquisition Unit, or 
prompts the operator for data regarding boiler voltage, steam 
temperature and pressure, coolant flow, and coolant 
differential temperature. 

From this data the heat transfer rate can be calculated. 



Q=mCp{T^-T^) 



(19) 



Then the overall heat transfer coefficient is calculated: 



Q 



AALMTD) 



( 20 ) 
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where : 



LMTD= 




( 21 ) 



Since the desired output is outside heat transfer coefficient, 
the principle of thermal resistances in series is used, where 
the tube wall thermal resistance is written as: 



DRPALL contains a computer code for the Modified Wilson 
Plot Technique to determine the inside and outside heat 
transfer coefficients. As described by Cobb [Ref. 8], the 
Modified Wilson Plot Technique uses the overall heat transfer 
coefficient to find the inside and outside heat transfer 
coefficients using assumed forms for them and following an 
iterative technique. Since the data were taken using the 
Petukhov-Popov correlation on the cooling water side[Ref. 10], 
the heat transfer coefficients were assumed to be: 




( 22 ) 



and the overall thermal resistance is given by: 




(23) 
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(24) 



*o = « 



\ifD,ATf 

k 



1/4 



*i=C, 



Q 



(25) 



where ; 

— RePr 

Q= 



(26) 



e=[1.82log(i?e) - 1 . 64 ] 


(27) 


ki=l+3 .46 


(28) 



and: 



i<2=11.7+1.8Pr‘^/3 (29) 

The values of a and C,- are calculated in the code. In 
addition, DRPALL contains corrections to take into account 
frictional heating of the coolant, as well as the fin effects 
of the two mounted ends of the test tube. More information for 
the Program DRPALL is given in Appendix A. 
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2 . 



HEATMEYER 



"HEATMEYER" is a computer code originally written by 
Cobb [Ref. 8] and called HEATCOBB. HEATMEYER is a slightly 
altered version of HEATCOBB in order to allow an interactive 
input of tube parameters. This program is written in FORTRAN 
and is a codified version of the Rose model [Ref. 4], with one 
very important difference. Cobb [Ref. 8] modified the Rose 
model to take into account the effects of fin efficiency. The 
same fin efficiency equation used by the Beatty and Katz model 
[Ref. 2 ], was applied. 

All numerical values of outside heat transfer coefficient 
and enhancement, presented in this paper, that are attributed 
to Rose (modified) are determined by using this program. More 
information for the program HEATMEYER is given in Appendix B. 

3 . Tsu j imor i 

In 1993, Tsujimori [Ref. 11], produced computer codes 
which calculate outside heat transfer coefficients and 
enhancements (for a given temperature difference) for the 
models of Nusselt, Beatty and Katz, Adamek and Webb, and Honda 
et al. All numerical values of outside heat transfer 
coefficient and enhancement presented in this thesis, which 
are attributed to Nusselt, Beatty and Katz, or Adamek and 
Webb, or Honda et al., were determined by use of Tsujimori's 
codes. More information regarding the Tsujimori programs is 
given in Appendix C. 
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V. RESULTS AND DISCUSSION 



A. GENERAL DISCUSSION 

Data were taken as described in Chapter IV, with two runs 
being done on each tube: one at atmospheric pressure, and 
another under vacuum conditions. Short form printouts of the 
data as taken and processed by program DRPALL are included in 
Appendix D. 

The names of the data files give information on the tube 
type and configuration, as well as the type of operation. The 
first two letters of the file name tell which type of tube 
material was used. For example, "ss” means stainless steel, 
and "cn" means copper-nickel. The numerical values in the file 
name represent the fin height of the tube where "15" means a 
fin height of 1.5mm, "125" means 1.25ram, "1" means 1mm etc.,. 
Finally, if the file name ends with an "A", that means the 
experimental data were taken at atmospheric pressure, vice a 
vacuum. Any file that ends with an "R" means that an original 
run had been terminated because of equipment problems , and 
that the run had been repeated. 

Any time experimental data are taken, experimental 
uncertainty becomes an important concern. Appendix E contains 
the program used to predict the uncertainty for any given run, 
as well as a brief explanation of the logic used. Appendix E 
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also contains the uncertainty analyses for all of the data 
runs . 

Related to uncertainty is the issue of repeatability. 
Consistency of experimental results is very important. In 
other words, it is vital that the data taken reflect the way 
tubes transfer heat, not the way the author collected his 
data. To demonstrate repeatability. Table III is a comparison 
of data taken by Cobb [Ref. 8] and the author for two tubes of 
identical dimensions ( 1mm fin height, 1mm fin thickness, and 
1.5mm fin spacing ) at vacuum. 

Another indication of repeatability is how the data from 
one tube compares with that of another, ie, are there any 
trends or does the data seem entirely random? As demonstrated 
in the plots to follow, there are some very clear trend which 
help establish the repeatability of any one individual data 
run. 



TABLE III. COMPARISON OF INDEPENDENT RUNS OF FINNED TUBES 



1 TUBE 

MATERIAL 


Ci 


alpha 


ENHANCEMENT I 
(delta T) 1 


copper- 








nickel 


2.33 


1.07 


1.32 


(Cobb) 
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copper- 

nickel 

(Meyer) 


2.68 


1.06 


1.30 


% difference 


13.1 


1.5 


1.5 


copper 


2.99 


1.50 


1.85 


(Cobb) 








copper 


2.87 


1.51 


1.86 


(Meyer) 








% difference 


3.9 


0.5 


0.5 



B. HEAT TRANSFER COEFFICIENT VS. TEMPERATURE DIFFERENCE 

Figures 5 through 12 are plots of the outside heat 
transfer coefficient versus film temperature difference where 
the temperature difference, again, is defined as the 
difference between the saturation temperature of the steam and 
the outside wall temperature of the test tube calculated at 
the base of the fin. Figure 5 also shows some sample 
uncertaint bars as determined in Appendix E. Two points 
immediately make themselves clear: 

1. Improvement of Enhanced Over Smooth Tube Performance 
For two tube materials, copper and aluminum, data were 
taken on smooth tubes with the same outside diameter as the 
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Ho v3. delta T for Aluminum tubes at vacuum 
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Figure 6 Experimental Results of Ho Vs. 

Temperature Difference for Aluminum 
Tubes at Vacuum 
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Temperature Difference for Copper-Nickel 
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diameter of the finned tubes at the base of the fins (ie the 
root diameter) . Exactly as one would expect, there is a marked 
increase in the heat transfer of the integral- 
fin tubes when compared to the smooth tubes. These effects can 
be seen tnp&cguoesCdpdvi&t^isiiAylOn Tube Performance 

When comparing the data for high conductivity 
materials, such as copper or aluminum, against the performance 
of low conductivity materials, such as copper-nickel or 
stainless steel, it becomes apparent that the conductivity of 
the material plays a large role in tube performance. There is 
a very definite trend established that as thermal conductivity 
decreases, so does heat transfer performance. The stainless 
steel plots in particular, (Figures 8 and 12) demonstrate that 
beyond fin heights of 0.5mm for vacuum, and 0.75mm for 
atmospheric, the effect of the low conductivity is so 
significant (ie, low fin efficiency) that the heat transfer 
coefficient does not increase with fin height. 

In fact, beyond these critical fin heights, the heat 
transfer coefficient decreases with fin height. This can be 
explained by the fact that, as described previously in Chapter 
I, as fin height increases, not only is fin efficiency 
reduced, but, the amount of tube that is flooded increases, 
reducing the amount of tube surface for effective condensation 
to occur, and therefore decreasing the outside heat transfer 
coefficient. 
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C. COMPARISON OF DATA WITH PREDICTIVE MODELS 



Figures 13 through 20 are plots of outside heat transfer 
coefficient against temperature difference for the 

experimental data and five predictive models. This is done 
for tubes of a fin height of 0.75mm. The models are those of 
Adamek and Webb [Ref. 3], Honda et al. [Ref. 5], Beatty and 
Katz [Ref. 2], modified Rose [Ref. 4], and Nusselt [Ref. 1] . 

The Nusselt model is for a smooth tube vice a finned tube 
and is only included to provide an indication of the 

enhancement achieved by using finned tubing. 

There are two models which seem to consistently predict 
tube performance reasonably well. They are the models of Rose 
(modified) [Ref. 4], and Beatty and Katz [Ref. 2]. 

The Beatty and Katz model, which, while reasonably 
accurate, consistently over-predicts the experimental 
performance of the integral-fin tubes. This is due to the fact 
that Beatty and Katz neglected the effects of surface tension. 
In fact, the Beatty and Katz model clearly is more accurate 
for the atmospheric runs than it is for the vacuum runs. This 
is because the atmospheric runs are conducted at 100 degrees 
C (vice 48.7 C during vacuum conditions) where the condensate 
surface tension is reduced. 

The modified Rose [Ref. 4] model appears to be overall the 
most accurate model, although it tends to under-predict the 
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experimental tube performance. Of course, accuracy coupled 
with conservatism can be a very desirable design 
characteristic. Actual values of enhancement as predicted by 
modified Rose, experimental enhancement, and the percent 
difference between the two, will be presented later in tabular 
form. 

As Figures 13 through 20 show, the Adamek and Webb [Ref. 
3] model tends to excessively over-predict the performance of 
integral-fin tubes. Though the model displays the correct 
trends, the relative inaccuracy and complexity compared to the 
modified Rose model, would tend to render the Adamek and Webb 
model unusable. 

The Honda et al. [Ref. 5] model demonstrates the ability 
to be extremely accurate, but its predictions vary widely as 
the model steps through its different sub-cases (the wide 
changes in outside heat transfer coefficient predicted by the 
Honda model do not seem to be borne out by the experimental 
results) . Again, the complexity and often inaccuracy of the 
Honda model makes other models such as modified Rose, more 
appealing. The inaccuracies of the Adamek and Webb [Ref. 3] 
and Honda et al. [Ref. 5] models may be due to errors in the 
codes established by Tsujimori [Ref. 11]. 

D. ENHANCEMENT VS. FIN HEIGHT 

Figures 21 and 22 are plots of the experimental 
enhancement ratio versus fin height for all four tube 
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Figure 22 Experimental Results of 
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materials. The enhancement is defined as the ratio of 
experimentally found outside heat transfer coefficient at a 
given temperature difference, over the outside heat transfer 
coefficient for the same temperature difference as predicted 
by Nusselt. There are three major points which can be derived 
from these plots: 

1. Smooth Tube Performance 

For the copper and aluminum smooth tubes (ie fin 
height equal to zero) , one can see a slight enhancement over 
that predicted by Nusselt. This is due to the fact that 
contrary to Nusselt 's assumption of a quiescent vapor, there 
is a downward vapor velocity associated with the experimental 
data (approximately 2 m/s for vacuum runs and 1 m/s for the 
atmospheric runs) . This vapor velocity tends creates a shear 
force that thins the condensate film and enhances heat 
transfer. 

2. Effect of Fin Height on Enhancement 

Again, particularly for high conductivity materials, 
as fin height increases, so does performance. For example, for 
copper and aluminum tubes, one can see an increasing 
enhancement up to a fin height of 1.5mm, and the data appear 
to demonstrate that a further increase in enhancement may 
occur if fin height is further increased. However, this is not 
so for low conductivity materials as discussed in the next 
section. 
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3. Effect of Conductivity on Enhancement 

Low thermal conductivity materials severely reduce 
enhancement. As can be seen in Figures 21 and 22, raising the 
fin height would not necessarily result in further enhanced 
performance. Even for a material with an intermediate thermal 
conductivity, such as copper-nickel (see Figure 21) , beyond a 
fin height of about 0.75mm, there is little increase in the 
enhancement. For stainless steel, the enhancement decreases 
for a fin height above 0.5 - 0.75mm, depending on the 
operating conditions. 

In the present study, the minimum fin height used was 
0.5mm. For stainless steel, it is observed that under vacuum 
conditions, the enhancement peaks at a fin height of 0.5mm and 
decreases for larger values. A recent work by Jaber and Webb 
[Ref. 6], shows that for titanium tubes, which have a 
conductivity near that of stainless steel, the enhancement 
increases with increasing fin height of 0.28 and 0.43mm. It 
appears that for such tubes, 0.5mm fin height would result in 
an optimum performance. However, more experimentation with 
lower fin heights is required before any firm conclusions can 
be reached. 

E. COMPARISON OF ENHANCEMENT WITH THE ROSE (MODIFIED) MODEL 

Figures 23 through 30 are plots of enhancement versus fin 
height and compare the experimental data to the predictive 
results of the modified Rose model. Note that for all the 
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Enhancement 



Enhancement vs. fin ht for Aluminum tubes at vacuum 




Figure 24 Experimental Results of 

Enhancement Vs. Fin Height for 
Aluminum Tubes at Vacuum 
with the Rose Model 
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Enhancement vs. fin ht for Cu-Ni tubes at vacuum 




Figure 25 Experimental Results of 

Enhancement Vs. Fin Height for 
Copper-Nickel Tubes at Vacuum 
with the Rose Model 
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Enhancement vs. fin ht for Stainless steel tubes at vacuum 




Fin ht (mm) 

Figure 26 Experimental Results of 

Enhancement Vs. Fin Height for 
Stainless Steel Tubes at Vacuum 
with the Rose Model 
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Enhancement vs. fin ht for Copper tubes at atmospheric pressure 
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Enhancement 



Enhancement vs. fin ht for Cu-Ni tubes (atmospheric) 




Figure 29 Experimental Results of 

Enhancement Vs. Fin Height for 
Copper-Nickel Tubes at Atmospheric 
Pressure with the Rose Model 
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plots, the modified Rose model demonstrates a reasonable to 
very good predictive capability. The only glaring shortcoming 
of the Rose model is its inability to predict the performance 
peaks of stainless steel at low fin heights (Figures 26 & 30) . 

Surprisingly, even though the original Rose model was 
developed using experimental data for copper tubes at 
atmospheric pressure, the modified Rose model works well for 
all tube materials. In addition, one might expect, that the 
modified Rose model would work best for copper tubes at 
atmospheric pressure, when in fact, this is not the case. This 
may be at least partially explained by recognizing that the B 
coefficients for the Rose model were determined without taking 
into account fin efficiency. Adding a fin efficiency to create 
the modified Rose model would then make the coefficients 
incorrect since they essentially include the effects of copper 
fin efficiency, assumed to be unity. Accuracy of the modified 
Rose model improves for conductivities less than that of 
copper, probably because the effects of fin efficiency become 
increasingly predominant. 

Table IV. compares enhancement for a given experimental 
data run to the average enhancement as predicted by Rose 
(modified) for the same film temperature difference. 

Note that with very few exceptions, the modified Rose 
model was able to predict the experimental data with good 
accuracy. The few exceptions may be more an indication of 
experimental error than of problems with Rose's (modified) 
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model. The potential of the modified Rose model warrants more 
experimental data to further establish its validity. 



TABLE IV. EXPERIMENTAL AND ROSE MODEL ENHANCEMENTS 



1 TUBE TYPE 


EXP 


ROSE 

(MODIFIED) 


% DIFF. 


CU5 


1.53 


1.41 


7.8 


CU75 


1.65 


1.60 


3.0 


CUl 


1.86 


1.82 


2.1 


CU125 


2.10 


2.04 


2.2. 


CUl 5 


2.16 


2.21 


2.3 


AL5 


1.27 


1.39 


9.4 


AL75 


1.65 


1.57 


4.8 


‘2.3 


1.73 


1.79 


3.§ 


AL125 


1.69 


1.96 


15.0 


AL15 


1.91 


2.09 


9.4 


CN5 


1.17 


1.S3 


3.0 


CN75 


1.28 


1.37 


7 . 0 


CNl 


1.30 


1.49 


14.6 


CN15 


1.38 


1.58 


14.5 
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TUBE TYPE 


EXP 


ROSE 

(MODIFIED) 


% DIFF. 


SS5 


1.20 


1.02 


15.0 


SS75 


1.12 


1.02 


8.9 


SSI 


0.96 


1.03 


6.7 


SS125 


0.91 


1.02 


12.1 


SS15 


0.96 


1.01 


5.2 


CU5A 


2.11 


1.47 


30.3 


CU75A 


2.14 


1.70 


20.5 


CU125A 


2.44 


2.16 


11.5 


CU15A 


2.60 


2.34 


10.0 


AL5A 


1.30 


1.44 


10.8 


AL75A 


1.75 


1.66 


5.1 


ALIA 


2.14 


1.91 


6.1 


AL125A 


1.96 


2.08 


6.1 


AL15A 


2.24 


2.23 


0.4 


CN5A 


1.32 


1.35 


2 . 3 


CN75A 


1.52 


1.48 


2.6 


CNIA 


1.61 


1.61 


0.0 
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TUBE TYPE 


EXP 


ROSE 

(MODIFIED) 


% DIFF. 


CN15A 


1.83 


1.72 


6.0 


SS5A 


1.36 


1.10 


19.1 


SS75A 


1.44 


1.13 


21.5 


SSIA 


1. 14 


1.15 


0.9 


SS125A 


1.17 


1.15 


1.7 


SS15A 


1.10 


1.13 


2.7 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Experimental data were obtained for steam condensation on 
integral-fin tubes made of copper, aluminum, 90/10 copper- 
nickel, and 316 stainless steel at both atmospheric and vacuum 
conditions. The tubes used had a root diameter of 13.88mm, a 
fin thickness of 1.0mm, a fin spacing of 1.5mm and fin heights 
ranging from 0.5mm to 1.5mm, in 0.25mm increments. From this 
data, the following conclusions can be made: 



1. Reliable, repeatable data have been obtained, on the 
performance of integral-fin tubes of varying materials 
and fin heights. 

2. For high conductivity materials, such as copper or 
aluminum, as fin height increases so does the 
enhancement of performance. 

3. For low conductivity materials, such as stainless 
steel, the effect of increasing surface area for heat 
transfer by raising fin height, is negated by both the 
poor fin efficiency, and the increased flooded area of 
the tube, resulting in a decrease in heat transfer 
performance. 

4. Of the examined predictive models, the modified Rose 
model seems to be the most accurate. This is despite 
the fact that his empirically determined coefficients 
were found only with data for a copper tube at 
atmospheric pressure. 
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B. RECOMMENDATIONS 



1. Use the results from tubes tested in this work and in 
future work, to evaluate the B coefficients in the 
modified Rose model to determine if the B values need 
to be changed. 

2. Test tubes at a fin height of 1mm with a fin spacing 
ranging from 0.5mm to 2.0mm to find a spacing which 
maximizes heat transfer enhancement for each tube 
material. 

3. Test tubes at a fin height of 1mm with a fin thickness 
ranging from 0.25mm to 1.5mm to find a fin thickness 
which maximizes heat transfer enhancement for each 
tube material. 

4. Using the results from 2 and 3, find the ideal fin 
configuration which maximizes heat transfer enhancement 
for each tube material. 

5. Experimentally determine how changing the root diameter 
of a tube changes the results in 4 . 

6. Continue with the computer upgrade in progress, to 
ensure faster, more timely analysis. 

7. Install a sight glass defogger on the test apparatus to 
enable the operator to easily visualize the tube during 
testing. 

8. Install a throttle valve to more precisely regulate 
the cooling water flow through the test tube. 
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APPENDIX A. 



PROGRAM DRPALL 



The computer program DRPALL, is a program written in HP 
Basic 3.0 which drives the HP 3497 Data Acquisition Unit. 
DRPALL takes the raw data, and using the Modified Wilson Plot 
Technique, calculates the test tube outside heat transfer 
coefficient. DRPALL also takes into account frictional heating 
of the test tube coolant, as well as tube end effects (ie it 
considers the fact that the two ends of the test tube act like 
fins) . 

More information on program DRPALL can be obtained by 
contacting : 



Prof. Paul J. Marto, Code ME/Mx 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey Ca. 93943-5002 



68 



APPENDIX B. 



PROGRAM HEATMEYER 



HEATMEYER is the program which predicts the outside heat 
transfer coefficient, and enhancement of integral-fin tubes 
based on the modified Rose model [Ref. 4]. HEATMEYER is a 
slight alteration of Cobb's HEATCOBB [Ref. 8]. More 
information on program HEATMEYER can be obtained by 
contacting: 

Prof. Paul J. Marto, Code ME/Mx 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey Ca. 93943-5002 
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APPENDIX C. 



TSUJIMORI COMPUTER CODES 



These codes, written by Tsujimori [Ref. 11] are written in 
the "C" computer language. There are a total of three 
individual programs. One program for Nusselt [Ref. 1] (as a 
reference), as well as Beatty and Katz [Ref. 2], one program 
for the Adamek and Webb [Ref. 3] model, and the last program 
for the Honda et al. [Ref. 5] model. 

All three programs are interactive and are written such 
that the user may specify the test tube parameters for any 
tube without having to alter the program. All three programs 
generate data files of heat transfer coefficient vs. 
temperature difference, as well as enhancement ratio vs. 
temperature difference or heat flux or fin spacing. For more 
information on the Tsujimori codes, contact: 

Prof. Paul J. Marto, Code ME/Mx 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, Ca 93943-5002 
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APPENDIX D. - EXPERIMENTAL DATA 
This Appendix has short form printouts, generated by 
program DRPALL, for all data runs taken. 
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APPENDIX E. 



UNCERTAINTY ANALYSIS 



When taking experimental measurements, error is always 
introduced. Though great care was used to ensure the accuracy 
of the data taken, there is no such thing as perfectly exact 
measurements. While the error introduced by any one particular 
measurement may be small, the cumulative error introduced by 
all the measurements may become quite large. 

Uncertainty is defined as the estimated difference between 
the actual measured value, and the calculated one. Kline and 
McClintock [Ref. 12] developed a method to determine the 
uncertainty of an experimentally derived value. This value V, 
which is a function of many measured quantities ie, V = 
V(x^,X 2 ,X 3 , . . .x^) , has an uncertainty given by the formula: 






f U,1 


2 




2 


y,l 


2 


[— C'f 










dxj \ 




”J) 



1/2 



(30) 



where: 

Uy = the uncertainty in the dependant variable 
x^,Xj,...x^ = the measured independent variables 
U^,U 2 ,...U^, = measured variable uncertainty 

Georgiadis [Ref. 13], gives a complete description of the 
uncertainty analysis used. 
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The uncertainty analysis program used is given in this 
Appendix along with examples, and was a revision of Cobb's 
[Ref. 8]. 
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At I 

V oCUUPi 



i rt . . 

V ucQ 



w ; 

Heat FIU’X 

TUu c~*r/\St a 1 thcPTASl CGHClUC . 
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